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Abstract 

We provide an action for gauge theories discretized on simpli- 
cial meshes, inspired by finite element methods. The action is dis- 
cretely gauge invariant and we give a proof of consistency. A discrete 
Noether's theorem that can be applied to our setting, is also proved. 



1 Introduction 

The Standard Model describes fundamental particles and interactions, except 
gravity, as a quantum field theory with gauge group U(l) x SU(2) x SU(3), 
see e.g. [33] [M]- Lattice Gauge Theory (LGT) [3E]EI[3m * s a computational 
approach with good agreement with experimental data. It has proved par- 
ticularly useful for the SU(3) sector, concerning quarks and gluons, where 
perturbative methods fail. 

In LGT, the underlying space-time is discretized by a cubical lattice. 
Scalar fields are assigned degrees of freedom on vertexes. Gauge fields are 
assigned degrees of freedom on edges, representing parallel transport between 
vertexes (Wilson lines). Curvature is recovered from holonomies around 
squares (plaquettes) in the lattice (Wilson loops). These quantities are com- 
bined to form a discrete action (Wilson action) that is gauge invariant, with 
respect to gauge transformations at vertexes. The construction of a discrete 
gauge invariant action is a fundamental ingredient in LGT, and is the only 
part of LGT we are concerned with in this paper. 

There is interest in constructing discrete gauge invariant actions on other 
underlying geometries, in particular simplicial decompositions of space-time. 
This has already received considerable attention^, starting with [3] [15] . For 
developments motivated by non- commutative differential geometry, see j!3j . 
The approach we propose in this paper results in an alternative prescription 
and is inspired by finite element methods, as they have been developed for 
Maxwell's equations, corresponding to the U(l) sector. 



1 We were blissfully ignorant about these, when submitting the first draft. 
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The most successful finite elements for Maxwell's equations are those in- 
troduced in [23], generalizing [2B]- As remarked in [3], lowest order Raviart- 
Thomas-Nedelec elements, correspond to Whitney forms [32] [35]. Based 
on this connection, a finite element exterior calculus has been developed 
[IE] [IS] [6] pQ . We refer to [2] for a recent review, relating the analysis de- 
veloped with engineering problems in mind, to discrete Hodge theory, as in 

Here, we tap into this large body of work, continuing our investigations 
of Lie algebra valued Whitney forms [10] on the one side and LGT [8] [7] on 
the other. The point of view we develop, uniting these strands, is that LGT 
should define an action on the space of Lie algebra valued Whitney forms, 
which is close to the restriction of the Yang- Mills action, yet invariant under 
some discrete gauge transformations. 

The Yang-Mills action is of course gauge invariant, but the space of Lie 
algebra valued Whitney forms is not, resulting in a violation of local charge 
conservation. The analogue for electromagnetics would be that the discretiza- 
tion would violate electric charge conservation. A discrete gauge invariance 
on the other hand, should give local charge conservation by discrete Noether's 
theorems, mimicking the continuum theory. A key ingredient in any conver- 
gence analysis is to estimate the error between the restriction of the Yang- 
Mills action and the new one. In the context of finite element methods, such 
consistency errors are analysed in the framework of variational crimes [3T] . 
see H2] §26 - 29. 

For gauge group U(l), a consistent and gauge invariant action on sim- 
plicial meshes can be obtained through a relatively simple construction and 
some applications can be completely analyzed [9]. A key ingredient is a notion 
of mass lumping [T7]. Because of limitations in the scope of mass lumping, 
we don't expect this method to extend to the full Yang-Mills action. 

By a more elaborate procedure, we define here, a consistent and gauge 
invariant discrete action for Yang-Mills theories on simplicial meshes. The 
simplexes are not congruent, so the metric must enter the formulas in a non 
trivial way, contrary to the cubical lattices of standard LGT. The metric 
on a simplex defines a mass matrix for the Whitney two-forms. Wilson 
loops associated with the faces of the simplexes are used to represent the 
curvature covariantly. Whereas standard LGT sums individual contributions 
from faces, we sum over simplexes, in which cross terms between the different 
faces appear. These couplings between Wilson loops are weighted by the mass 
matrix coefficients and made gauge invariant by discrete parallel transport 
between origins. 

One advantage of the proposed formulas is to allow local mesh refine- 
ments, useful for instance to efficiently represent singular fields. The con- 
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sistency proof we provide covers such meshes, indicating robustness with 
respect to mesh geometry. Another advantage is to accommodate variable 
metrics as defined by Regge calculus [29] , see Remark |3j 

This paper serves to introduce the formalism, give the definition, prove 
consistency and propose a discrete Noether's theorem. Preliminary applica- 
tions to quantum field theory, including numerical results, have been reported 
elsewhere [IB] , while this paper was under review. It is organized as follows. 
Section [2] contains definitions pertaining to connections and curvature, as 
well as Whitney forms. Section [3] includes the definition of the proposed 
discrete Yang-Mills Lagrangian and some comments. Section H] contains the 
proof of consistency in a finite element sense. Finally section [5] contains the 
discrete Noether's theorem we introduce. 

2 Prerequisites 

Yang-Mills action A standard reference for connections and curvature is 
|2Uj . Here, we use notations as in [TO], which also contains a more com- 
prehensive presentation of Lie algebra valued differential forms. The main 
ingredients are as follows. 

Choose a compact Lie group G with associated Lie algebra q. For simplic- 
ity we suppose that G is a subgroup of the complex unitary n x n matrices, 
for some n. Typically an element of G will be denoted G and an element of 
q will be denoted g. The Hermitian conjugate of a matrix g is denoted g H 
and the real scalar product is: 



When no confusion is possible with scalars, the unit matrix is denoted 1 and 
the zero matrix is denoted 0. 

Let S be a bounded domain in m-dimensional Euclidean space (in ap- 
plications m = 2,3,4 would be the most common). The space of smooth 
k- forms on S is denoted Q k (S). The space fl k (S) ® Q can be identified with 
the space of smooth q- valued fc-forms on 5*. The bracket of Lie algebra valued 
forms is determined by: 



where u, u' are real valued differential forms and g, g' are elements of q. 

A smooth connection one-form on S is an element A £ £g> q also 

called gauge field. Its curvature is J 7 (A) £ Q 2 (S) ® q defined by: 




(1) 




(2) 



F(A) = dA + l/2[A,A]. 



(3) 



3 



We will use such forms with less regularity, typically in some Sobolev space. 

Gauge transformations of connection one-forms are associated with func- 
tions Q : S — > G, and defined by: 

Q Q {A) = QAQ- 1 - (DQ)Q- 1 . (4) 

One has: 

HG Q {A))=QHA)Q~ l . (5) 
The Yang-Mills action is given by: 

S(A) = I \F{A)\ 2 . (6) 
Js 

Since the adjoint representation is unitary, this action is invariant under 
gauge transformations. 

Whitney forms We refer to [H] pQ for surveys on Whitney forms [32] [35] 
in a finite element guise. For aspects relating to differential geometry and 
algebraic topology, one can consult [27]. The following is a summary of 
results needed and serves mainly to introduce notations. 

Let T be a simplicial complex, spanning the domain S. The set of k- 
dimensional simplexes in T is denoted T k . Simplexes of dimension 0, 1 and 
2 are referred to as vertexes, edges, faces respectively. Generic labels for 
edges and faces will be e and / respectively. The symbol T can be used for 
simplexes of any dimension. In the presence of several vertexes we denote 
them by i, j, k, I. We suppose an orientation has been chosen for each simplex 
in T. 

Let W k (T) be the space of Whitney A;- forms on T. We also denote by 
W k (T) the space of Whitney fc-forms on a simplex T. The canonical basis of 
W k (l~) is denoted (At), T ranging over the set T k of k- dimensional simplexes 
in T. Explicitly, the 0-forms are spanned by the barycentric coordinate maps. 
For any vertex i G T°, \ is the piecewise affine map taking the value 1 at 
vertex i and at the other vertexes. For a k > 1 and a fc-dimensional simplex 
T G T with vertexes io, ■ ■ ■ ,ik, ordered according to the chosen orientation 
of T, we have: 

k 

K-io = A T = A;! ^(-1^'A^dA^ A ... dA*. ... A dA lfc . (7) 

j=0 

The hat signifies omission of this term. We will only use 0-, 1- and 2- forms. 
When i, j are vertexes of an edge, the associated Whitney 1-form is: 

Xji = AjdAy — AjdAj. (8) 
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When k are vertexes of a face, the associated Whitney 2-form is: 

X kji = 2(X i d\ j A dAfc - XjdXi A dA fc + X k dXi A dA^). (9) 

The space of (real) /c-cochains consists of the functions that assign a real 
number to each fc-simplex: 

C k (T)=R Tk . (10) 

The coboundary operator 5 : C k (T) —> C k+l (T) is defined as follows. The 
relative orientation o(T, T") between a simplex T £ T k + 1 an d a simplex 
rpi ^ -y-fc j g g j£ j g n0 ^. j n ^ e k ounc iary Q f anc [ when T" is in the 

boundary of T, the sign depending on whether it is outward oriented or not. 
For u £ C k (T), 5u £ C k+1 (T) is defined on any T £ T k+1 by: 



Nr= ^ o(T,T')u 2 v. 

One has 55 = 0. The de Rham map R k is defined by: 

Rk . | fi fc (S) C fc (T), 

M ^ (f T u) TeT k. 



(12) 



In this formula the /c-form w is (pulled back and) integrated on the k- 
simplexes T, taking into account orientations. By Stokes' theorem we have 
commuting diagrams: 

ft fc (s) — n k+1 (S) 



(13) 



R k+1 



C fc (T) — C k+1 (T) 
Since for T, T' £ T fe one has: 

At = Stt' (Kronecker delta), 



L 



T' 

the de Rham map induces isomorphisms: 

R k : W k {T) -> C fc (T), 

whose inverses are: 

C fc (T) -> W fc (T), 



(14) 



(15) 



(16) 
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Given u G W k (T) we denote by u. = R k u its associated cochain. Another 
useful notation is the following. If u G C k (T), one defines for any simplex 
TeT with vertexes io, ■ ■ ■ , ik '■ 

U i k -io = ±U T, (17) 

the sign depending on whether the ordering of the vertexes agrees with the 
orientation of T or not. For instance we can write if u G C 2 (T): 

(5u) k ji = u ik + u k j + Uji. (18) 

Let I k denote the interpolation operator onto Whitney /c-forms - it is the 
projection onto W k (T) determined by the identity R k I k = R k . Equivalently, 
for u G Q k (S) one has: 

I k u= J2(J T ^T£W k (T). (19) 

TeT k 

Interpolation commutes with the exterior derivative. 

Whitney forms are not smooth, but have enough regularity for the exterior 
derivative, in the sense of distributions/currents of Schwartz and de Rham, 
to be given by the simplex-wise definition (there are no Dirac measures on 
interfaces) . 



3 Definition 

We first make some remarks on Lie algebra valued Whitney forms and their 
associated cochains. The de Rham isomorhisms ffT5l) give isomorphisms: 

R k : W k {T) ® -> C k {T) ® S, (20) 

An assignment of a Lie algebra element to each fc-simplex will be called a 
Lie algebra /c-cochain. For k — 1 this goes as follows. Pick A G W l (T) ® Q- 
Attached to an edge with vertexes i, j and oriented from i to j, one has an 
element Aji G g obtained by writing: 

A = ^ AjiXji where A j{ = A (21) 

Thus ji denotes an edge in T 1 , oriented from the vertex i to the vertex j. 
With these notations, the cochain associated with A is: 

A. = (Aji)j i( z T i G C l {T) <8> Q- (22) 
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The pullback of the 1-form A G W l (T) ® Q to the edge jz is a constant 
1-form. Therefore parallel transport from i to j is given simply by: 

Uji = expi-Aji). (23) 

We suppose Uji to be close enough to 1 for the logarithm to be unambiguous. 
Then one is free to think in terms of Lie group elements Uji (close to 1) or 
Lie algebra elements Aji (close to 0). We use the conventions: 

U ij = Ujl 1 and Uu = l, (24) 

which correspond to: 

Aij = -Aji and A i{ = 0. (25) 

A discrete gauge transformation is associated with a choice of Lie group 
elements G G, one for each vertex i G T°. One then transforms A G 
W 1 (T) ® by operating on its parallel transports. Namely, if Uji are the 
parallel transports of A, the parallel transports of its image A' G W 1 (T) <8> 
will be U^ defined by: 

= GjUjiG' 1 . (26) 

With this choice of gauge transformations, we will construct a gauge invariant 
approximation of the "true" Yang-Mills action on simplexes T G T m of 
maximal dimension, which we recall to be defined by: 

St{A) = J t \HA)\ 2 - (27) 

In our setting T inherits the Euclidean metric of the ambient space, but as 
already indicated one could use a Regge metric instead. The metric enables 
integration of scalar functions on T. It also gives, at each point x of T, 
a scalar product on alternating forms above x. The associated norm was 
denoted | - | in (EZD- 

Let M be the matrix of the L 2 (T) product on W 2 (T), equipped with the 
standard basis defined by ()9]). This matrix is indexed by the two-dimensional 
faces of T (and depends on their orientations). Explicitly, for two faces fo 
and fi in a simplex T, we put: 

M foh (T) = J \ h ■ X h . (28) 

where the scalar product of alternating forms is denoted (•). Notice that this 
matrix is not diagonal, which can be interpreted as an interaction between 
neighboring faces. 
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Given parallel transports U G C 1 (T, G), the discrete curvature associated 
with a face / G T 2 with vertexes fc is defined in analogy with square 
Wilson loops [36] by: 

Fkji = UikUkjUji. (29) 

In other words, one considers the holonomy around the boundary of the face 
/ of the 1-form A G W 1 (T) <8> Q related to U by (I231) . This formula depends 
on the ordering of the vertexes and locates the curvature at the vertex i. The 
curvature at vertex j is obtained by permuting indices and satisfies: 

Fikj = UjiF^jiUij. (30) 

This gives a formula for parallel transport of curvature from % to j. Con- 
cerning orientation of a given face, we also notice that the definition (129~|) 
implies: 

F jki = F kji- ( 31 ) 

Under gauge transformations this curvature behaves as follows. If F G 
C 2 (T, G) is the curvature associated with holonomies U G C : (T, G) by (1291 . 
and U is transformed by G G C°(T, G) to U', according to (J26]) , then the 
curvatures F' of U' are: 

F'kji = GiF^iG' 1 . (32) 

When / is a face with vertexes i, j, k, which is ordered as i — > j — > k and we 
choose to locate the curvature at i, we put: 

Ff = F kji- (33) 

This formula defines the curvature of a pointed oriented face f. For a pointed 
face /, its distinguished point is denoted / and called its origin. 

We now propose the following definition of a discrete action for lattice 
gauge theory on simplexes, as an alternative to (1271) : 

Definition 1. We define: 

S' T (A) = ^M /0/l (T) mttr (U fifo (l - F%)U hh {l - F h )). (34) 

/o/i 

In this formula we sum over pairs of faces /o, f\ of T, each one having an 
orientation and a distinguished point. We have incorporated the parallel 
transport determined by A, between the distinguished points fo and f\. 

We can state: 

Theorem 1. The action S' T is discretely gauge invariant, with respect to 
transformations of the form /[2b]) . 



S 



Global actions are obtained by summing the contributions of each maxi- 
mal simplex in T m , so that for A G W 1 (T) <g> £j: 

S'(A) = £ S 't(A). (35) 

Remark 1. Define, for any U G C 1 (T, G), subject to ff24|) and any F G 
C 2 (T, G) (not necessarily the curvature of £/!): 

£ T (C/,F) = ^M /()/l (T)^etr (f/ /i/o (l - - F h )). (36) 

/o/i 

For fixed U, £t{U,F) can be interpreted as an L 2 (T) norm squared of F. 
This norm depends on U, which contrasts with the fact that the usual L 2 (T) 
norm, as it appears in (|2"7|) is independent of the gauge. 

However Ct(U,F) is invariant under transformations (U,F) t— > (U',F'), 
associated with some G G C°(T, G) by (I2S|) and (I32"j) . 

Remark 2 . It is most natural to compute a norm in the Lie algebra. Thus in 
the definition of the discrete actions, the terms of the form [Ff — 1) should 
be considered as approximations of log Ff that are more readily computable. 
For a general Lie group one could use: 

S' T {A) = J2M fofl (T)Ad(U Afo )\og(F fo ).\og(F fl ). (37) 
foh 

Here Ad : G — > End($j) is the adjoint representation. Recall that for a given 
U G G, Ad(U) : S — ^ is the tangent map at unity, of the automorphism of 
G mapping an element G to UGU^ 1 . The scalar product on g, denoted here 
with (•), should make the adjoint representation unitary. 

Remark 3. The proposed method can be combined seamlessly with Regge 
calculus [29] (see [5] [11] for finite element interpretations). In Regge calculus 
the metric of a given simplex is determined by the edge lengths and yields 
a mass matrix, as in the adopted setting. It seems that only minor modifi- 
cations are necessary for the consistency proof we will provide, to cover the 
case where the local metrics are Regge metrics interpolating a smooth one. 
Notice that this action uses just edge lengths for the metric and values of 
the fields at vertexes or edges, never vertex coordinates, so that the method 
is "coordinate free". 

Scalar fields. We include a definition of a discrete action for certain so- 
called scalar fields. We will not prove consistency for it here. 
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Let V be an inner product space on which G acts unitarily. The action 
is denoted simply (G, v) !->■ Gv. Likewise the associated action of g on V is 
denoted (g,v) H- gv. 

Let V denote the canonical flat connection acting on sections $ : S — > V. 
Given A, the action to approximate is: 

S T (A,$) = [ |V$ + A$| 2 . (38) 
Jt 

Let then $ G M^°(T) ® V be a discrete scalar field. We can write: 

$ = $^ with = $(z). (39) 

Concerning the cochains associated with $ and V$ G W /1 (T) ® V we have, 
by (USD: 

V$= <I>,)A„. (40) 

The mass matrix for Whitney 1-forms is also denoted by M and is indexed 
by oriented edges. Thus: 

M eoei (T)= f A eo -A ei . (41) 
Jt 

For an oriented edge e we denote its origin by e and its target by e. 
As a discrete action we propose to use: 

S' T (A,$) = ^M eoei (T)C/e lgo ($e - £/ eo $e ) • - U^), (42) 

eoei 

where the scalar product is that of V. 

Recall that under discrete gauge transformations associated with G\ G G, 
the parallel transports U transform by ( |26|) . The corresponding transforma- 
tion of $ is, at the level of cochains: 

-)• Gi^i. (43) 
It is readily checked that S'(A, $) is discretely gauge invariant. 

Conventional LGT For comparison, we recall the usual definition of LGT 
on cubical meshes. One attaches a discrete parallel transport Uji to any two 
vertexes i, j of the grid linked by an edge, with the preceding constraint 
Uji = U~, . A face / of this mesh is then a square with four vertexes, called 
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plaquette. Given a choice of orientation and origin, these four vertexes can 
be labelled /o,/i,/2j/3- Barring the coupling constant, the action is then 
defined by a sum over faces: 

£<Ketr(l - U fof3 U hh U hfl U flfo ), (44) 
/ 

Remark that the action is independent of the choice of origin and orientation 
of the faces. 

Thus in standard LGT, one sums over faces, and the contribution of 
each face is discretely gauge invariant, under transformations ( 126ft . On the 
other hand, for simplicial meshes, we propose to sum over maximal simplexes 
(tetrahedrons in dimension 3), with gauge invariant cross terms between 
neighboring faces. The counterpart for cubical meshes would be to sum over 
cubes, inside which plaquettes are coupled two by two. From this point of 
view, standard LGT uses just diagonal terms. This is also why the scalar 
product (JT]) is more apparent in (134)) than in (1441) . One can compare with the 
fact that the Yee scheme [H7] can be deduced from a finite element scheme 
via mass lumping [21] (whereby mass matrices are approximated by diagonal 
matrices in a consistent way). 

That the continuum limit of the LGT action is the Yang Mills action, is 
usually argued on the basis of Taylor expansions and a couple of terms in the 
BCH formula, e.g. [5T] p. 786. To a numerical analyst these arguments would 
prove consistency in the finite difference sense. Consistency in the finite 
element sense is related of course, but not identical, putting emphasis on the 
choice of function theoretic norms, typically Sobolev norms. Transposing 
the finite element arguments we will give here for the simplicial case, to 
the cubical case (with tensor-product Whitney forms) would give a novel 
proof on the convergence of the Wilson action to the Yang-Mills action, in 
the continuum limit. We don't expect consistency proofs based on Taylor 
expansion techniques to carry over directly to the simplicial setting, since it 
is difficult for them to take into account mesh geometry. 

4 Consistency 

In this section we want to study the error committed, when approximating 
(T27|) by (l34|) . For this purpose we introduce two more discrete actions defined 
for A e W l (T) Cgijj, for simplexes T e T m . These discrete actions serve only 
to provide intermediate steps between St (A) and S T (A), aiming at clarifying 
the analysis. 
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First we define: 

St(A) = J t \I 2 HA)\\ (45) 

and remark that, by (flj). (fl9"j) and (j2"8"j): 

SUA) = £ M /0/l (D <Ke tr ( f fo J-(A) U j h J 7 (A)) . (46) 
foh 

The sum extends over pairs of faces fo, f\ of T. Second, given also a choice 
of origins of the faces, we define: 

SUA) = M foh(T) <Ketr ((1 - F /o ) H (l - F fl )). (47) 

foh 

We now evaluate the error in each of the approximations: 

S T ->■ S T -> S£ ->■ 5^. (48) 

The reasoning is partly similar to [9] , but we need to extend to non-commutative 
gauge group and abandon the convenience of mass-lumping. For definiteness 
we restrict attention to dimension m = 3 for the ambient space S. Maximal 
simplexes are then tetrahedrons. 

We suppose that we have a regular sequence of simplicial meshes T n of 
the domain S. The diameter of a simplex T is denoted Iit, and the biggest 
hr when T is in T n is denoted h n . We suppose that the sequence (/i n )neN 
converges to 0. Let denote the interpolant onto Whitney /c-forms associ- 
ated with the mesh l~ n . Let X n denote the space W 1 (T n ) <8> 0- For ease of 
notation put also I n = l\ and J n = I„. 

The following definition is the natural extension of [12] §28 to a non-linear 
setting. 

Definition 2. Let || ■ || denote the norm of some Banach space in which 
1 (S') (&q is dense, and containing each X n . We say that two actions S n and 
S' n defined on X n are consistent with each other, with respect to || • ||, if for 
all smooth A we have: 

sup \BS n (I n A)A' -DS' n (InA)A'\/\\A'\\ ^0 when n -> oo. (49) 

A'ex n 

More precisely if the above expression is 0(e n ), for some sequence e = (e„) — > 
0, we speak of consistency of order e. 
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If there is a constant C > (which may depend on A and the sequence 
(T n ) but not on n) such that quantities a n and b n satisfy a n < Cb n for all n, 
we write a n ^ b n or a n = 0{b n ). 

Consider a simplex T of dimension <i and let $ : T — > T be a scaling map 
of the form $(x) = /it^ + £/• We have, for any u e Q k (T): 

\\u\\lp(t) = h^ k+d/p \\^*u\\ hP( fy (50) 

Arguments based on this identity will be referred to as scaling arguments. 

For instance if we have a sequence of elements u n G X n which is bounded 
in L A (S) and e n is an edge in T n we deduce by scaling that we have bounds: 



u n \ 



< hT- (51) 



which gives: 

Kl^oc = max|<| r< hl/\ (52) 

This is enough to guarantee that the logarithm is unambiguous as required 
initially. 

On cochains, we consider norms of the following form (with no coeffi- 
cients): 

for u e C k {T) \u\ 2 = M 2 - (53) 

TeT k 

We will frequently use that on reference simplexes T, the cochain norm is 
equivalent to functional norms, as they appear in the right hand side of (1501) . 

We let A e <g> q be smooth. We put A n = I n A. Remark that for 

edges e in 7^ we have: 

A n e = (I n A) e = [a = 0(h e ), (54) 

J e 

and for faces / in T n we have: 



(5A:) f = (dI n A) f = fdA = 0{h}). 



(55) 



Step one We compare S and S 1 . We first remark: 
Lemma 1. We have: 

\\HA n ) - JnHA n )\\ mT) = 0{h T \\HA n )\\ mT) \ (56) 

and: 

\\JnHA n )\\v(T) ± \\HA n )\\v(T)- (57) 
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Proof. By scaling, knowing that T(A n ) lives in the space of Whitney forms 
of maximal polynomial order 2 ( [10] §3.3). □ 

Recall the formula: 

D\ A T(A)A' = dA' + [A, A']. (58) 

Since Whitney forms are stable under the exterior derivative, we have: 

DT(A n )A' - J n DT(A n )A' = [A n , A'} - J n [A n , A']. (59) 

Lemma 2. We have: 

\\[A n ,A'} - J n [A n ,A']\\ L 2 {T) = 0(h T \\A'\\ LHT) ). (60) 

Proof. Remark that the interpolation is exact on a tetrahedron if A n is con- 
stant on it. On a reference tetrahedron we can therefore write: 

||[y4 n ,v4'] — J n [A n , A']\\ L 2(f^ ^ IIV^IIlcx^IIA'Hj^j.). (61) 

The estimate on T then follows by scaling. □ 

Proposition 1. The actions S and S l are consistent of order h for the L 2 
norm. 

Proof. We have: 

DS T (A n )A' = [ T(A) ■ DT{A)A', (62) 
Jt 

and: 

DS^(A n )A' = [ J n T(A) ■ D J n T{A)A'. (63) 
Jt 

With a more compact notation we can evaluate the difference: 

| /( T - JT) ■ DTA' + / JT- (DJ-- JDT)A'\ ± h T \\T\\ L 2 {T) \\A'\\ L 2 {T) . 

1 (64) 

Summing these estimates for all the tetrahedrons and applying the Cauchy- 
Schwarz inequality gives the result. □ 
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Step two We compare now S 1 and S 2 . 

Lemma 3. Let a face f have vertexes k. We have: 

J^(A n ) = jdA n + ^[A n ,A n ], 

- A n -i- A n 4- A n 4- ~(\A n A n 1 -i- \ A n 4 n l -I- \ A n A n ~\} 

Proof. Remark that: 

JdA n = A- k + Al 3 + A^ 
Using formulas of the type: 

j Xji A \ k j = 1/6, 

one gets the second term. 
Proposition 2. We have: 

1 - Ff(A n ) = J T(A n ) + 0{h)). 

Proof. Ff(A n ) can be estimated with the help of the BCH formula and 
pared with the formula previously obtained for the right hand side. 

Using the same arguments as in the proof of Lemma [31 we get: 

Lemma 4. We have: 

D\ A (f f JT(A))A'= J f dA' + [A, A'} 

-A' ik + l/6([A%A' ik ]-[A n kj ,A' ik ])+ 
A'kj + l/6([-4"fe) A' k j] — [A^, A' kj ])+ 

4, + i/6(K J ,4,]-K fc ,4J). 

Proposition 3. We have: 

- D\ A F f (A)A' = D\ A j f J z (A)A' + 0(h f \6A' m \ + h}\A'.\). 
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Proof. Define an entire function 0, by setting, for z ^ 0: 

(j>(z) = (l-e-")/z. (75) 

Recall that: 

F kji = exp(-A ik ) exp(-A kj ) exp(-A ji ), (76) 

so that: 

-DF kji (A)A' = exp(-A ik )(j)(ad(-A ik ))A' ik exp(-A kj ) exp(-A^)+ (77) 
exp{-A ik ) exp(-A kj )</>(ad{-A kj ))A' kj exp(-A ji )+ (78) 
exp(-Afc) exp(-A kj ) exp(-A ii )0(ad(-A ji ))^ i . (79) 

Expand using <f)'(0) = —1/2 and rearrange to obtain, up to the announced 
error term, the previously computed right hand side. □ 

Proposition 4. The discrete actions S 1 and S 2 are consistent of order h 
with respect to the norm defined by: 



= \\dA\\ L * + \\A\\v. (80) 
Proof. Let fo and f\ be faces of a tetrahedron T in l~ n - Define: 

L n A> = D\ A=A «(f fo T(Arf h T(A))A>. (81) 

and: 

L' n A' = D\ A=An ((l - F /() (A)) H (1 - F fl (A)))A'. (82) 
Combining Propositions [2] and [3] gives: 

L' n A' = L n A' + 0(h 3 T \5A'.\ + 4|A'J). (83) 

We have, by scaling: 

\M Ml \^h^. (84) 
Insert it in the error term of (1831 and write: 

h 2 T \8A' m \ T + h%\A' m \ T =< 4 +1/2 ||dA'|| L2(r) + h 3 f 1/2 \\A'\\ L 2 {T) . (85) 

We sum over all tetrahedrons and apply a Cauchy-Schwarz inequality, 
remarking that: 

( h T) 1/2 1 h n- (86) 

TGT 3 

This concludes the proof. □ 
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Step three We compare S 2 and S'. 
Proposition 5. We have: 

UiiFkjiUu = Fkji + 0(11^), (87) 

and also: 

D\ A U H (A)F kji (A)Ua(A)A' = D\ A F kji (A)A' + 0{h T \5A',\ + h 2 T \A'.\). (88) 
Proof. For the first assertion we write: 

UiiFkjiUu — F^i = Uu(Fkji — l)Uu — (F^i — 1), (89) 
and conclude using: 

F^ ~ 1 = 0{h%). (90) 
For the second one we compute: 

D\ A U u F kji UuA' = U u V\ A F kji {A)AU a - U u [^d(-A U ))A' U , F kji ]U«. (91) 

On the right hand side, in the second term, we can replace F k ji by F k ji — 1 
and use again (I9"U|) . From this the second assertion follows. □ 

As in the previous paragraph we can deduce: 

Proposition 6. The discrete actions S 2 and S' are consistent of order h 
with respect to the norm defined by: 

\\A\\ = ||dA|| L2 + ||A|| L2 . (92) 

Conclusion Adding the three estimates proved in Propositions [TJ |4] and El 

we get: 

Theorem 2. The discrete actions S and S' are consistent of order h with 
respect to the norm defined by: 

\\A\\ = \\dA\\ L 2 + ||A|| L 2. (93) 

The arguments introduced also immediately show that, concerning the 
action itself, we have consistency of order h 2 . That is, if A is a smooth gauge 
potential, we have: 

S n {I n A)-S' n {I n A) = 0{h 2 n ). (94) 
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5 A discrete Noether's theorem 



In this section we propose an analogue of Noether's first theorem [26J, ex- 
pressed for discretizations over simplicial complexes, when the group acting 
on the fields preserves fibers, as defined below. This discrete result provides a 
discrete conservation law associated with discrete gauge invariance. While it 
does not capture the full power of the continuous one, it is sufficient to prove 
constraint preservation for evolution problems as in [8]. Recall that non- 
invariance of Lie algebra valued Whitney forms under discrete gauge trans- 
formations makes the standard action (|2"T|) problematic for the simulation of 
evolution, because constraints are not preserved [TU]. For electromagnetics, 
the conservation law associated in the continuum with gauge invariance, is 
nothing but electric charge conservation, of great physical significance. 

Related discrete Noether's theorems have been discussed in particular in 
[23] [22]. 

We suppose we have a finite simplicial complex T. The maximal dimen- 
sion of the simplexes in T is m. We write S <T to say that 5* is a subsimplex 
of T. If S is a simplex and i a vertex not in S, S + i is the simplex obtained 
by adjoining the vertex i to S. Conversely, if S is a simplex and i a vertex 
of S, S — i is the face of 5" opposite i. 

We suppose that on each maximal simplex T 6 T m we have attached 
fields $t of the form: 

$ T = ($ T (S))<j <T e 1[Vt(S). (95) 

That is, $t attaches a value in some space Vr(S) to each subsimplex S of 
T. We call Vr(S) the fiber above 5*. 

We suppose that we have a Lagrangian Ct attached to T, which is a 
function: 

C T : Y[V T (S) -»R. (96) 

S<fT 

In the following we fix a simplex T G T m . We suppose we have a one 
parameter group action which acts separately on each fiber Vr(S'): 

A T (S) : R -»■ Aat(V T (S)), (97) 

and for tel: 

A T [t]$ T = (A T (S)[t]$ T (S)) s<T . (98) 
We suppose that this group action leaves Ct invariant: 

Vt G R £ T (A r [t]$T) = C t {<$>t)- (99) 
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We define the (local) infinitesimal generators: 

fr(S)=d\ 1 ^ A T (S)[t]Q T (S). (100) 
and the (local) Euler-Lagrange functions: 

Er(S) = d\ s C T ($T), (101) 

and put: 

F T (S) = Et(S)Zt(S). (102) 

For each simplex S < T and each i G T \ S choose a number px(i, S) 
subject to the condition that, for any simplex S' of dimension at least 1: 

J> T (i,S"-i) = l. (103) 

ieS' 

Proposition 7. Define, for any vertex i G T: 

W T (i)=F T (i)+ J2 MhS)F T (S + i), (104) 

and for any two distinct vertexes i,j G T: 

V T {i,j)=F T {i)-F T (j)+ PT(i,S)F T (S+i)-p T (j,S)F T (S+j). (105) 

Then we have: 

(m + l)W T ({) = J2Mhj)- (106) 

Proof. In this proof, in which T is fixed, we drop the index T. The summation 
variables S 1 , 5" are subsimplexes of T. First we remark: 

E( F W+ E pM)F(S + z)) (107) 



mW(i) - E E P(hS)F(S + i). (10* 



then we remark: 



E(^(J')+ E PU,S)F(S + j)+ J2p(^ S ) F (S + i)) (109) 

E pu^-ms^i (no) 

En5")-^«- (HI) 
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From invariance of the Lagrangian we get : 

^F(50=O, (112) 

S' 

and this concludes the proof. □ 

In the applications we have in mind, if a simplex S G T is included in 
two maximal simplexes T,T' G T m we have Vr(S) = Vt'(S), and the global 
variable <J> has the property <&t(S) — &t'(S). When this happens for all 
choices S,T,T' such that S <T,T' G T m , we have a well defined fiber Vs 
above each S G T and the action S will be of the form: 

s = £t '][ v ^ ( 113 ) 

TgT m 5eT 

Moreover we suppose that the group action A acts separately on the fibers 
V(S), independently of any embedding into a maximal simplex T. In this 
setting we define the (global) infinitesimal generators: 

aS)=d\ t=0 A(S)[tMS), (114) 

the (global) Euler Lagrange functions: 

E(S) = d\ s jr($) = Et W> ( 115 ) 

T£T m :S<T 

and put: 

F{S) = E{S)£{S). (H6) 

We suppose finally that we have chosen the numbers pr{i, S) independently 
of T containing % and S. When i and S are not included in any simplex of T 
we set p(i, S) = 0. The preceding Proposition gives, by adding contributions 
from all maximal simplexes T: 

Proposition 8. Define, for any vertex i G T: 

W{i) = F{i)+ P(hS)F(S + i), (117) 

and for any two distinct vertexes i,jeT linked by an edge: 

v(i,j)= Yl Mhi)- (us) 

TeT m :i,jeT 

Then we have: 

(m + l)W(i)= V ^ti- ( 119 ) 

j-.i+jeT 1 
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In brief, equation (11191) expresses a weighted sum of (global) Euler-Lagrange 
functions applied to infinitesimal generators, as a discrete divergence. In- 
deed it is natural to think of V(i,j) as degrees of freedom of a vectorfield V. 
Choose any cellular complex dual to T, so that, in particular, the domain is 
covered by cells dual to the vertexes % G T°. Then V(i,j) is the flux from 
the cell dual to i into the cell dual to j, through the dual face of the edge 
ij G T 1 . The right hand side of (11191) is then the total flux leaving the cell 
dual to i, which is the natural degree of freedom for the divergence of V. The 
essential antisymmetry property V(i,j) = —V(j,i) guarantees that summing 
the discrete divergence over a union of top-dimensional dual cells, leaves only 
a boundary term. 

These considerations apply directly to the proposed simplicial gauge the- 
ory, for which moreover we have variables attached only to 0- and 1- sim- 
plexes. 
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